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We consider a model inspired by a metal break-junction hypothetically caught at its breaking 
point, where the non-adiabatic center-of-mass motion of the bridging atom can be treated as a two- 
level system. By means of Numerical Renormalization Group (NRG) we calculate the influence of 
the two level system on the ballistic conductance across the bridge atom. The results are shown to be 
fully consistent with a conformal fleld theory treatment. We find that the conductance, calculated 
by coupling Fermi liquid theory to our NRG is always finite and fractional at zero temperature, but 
drops quite fast as the temperature increases. 
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I. INTRODUCTION 

A number of transport measurements on organic and 
inorganic molecules bridged bewteen metallic leads have 
recently succeeded in revealing signatures of the molec- 
ular vibrational and motional degrees of freedom in the 
inelastic tunneling spectrum, and raised interesting the- 
oretical issues. Most notably, since in these nanosized 
devices the time scales of the nuclear dynamics may 
be comparable to those involved in the electron tunnel- 
ing, non-adiabatic quantum effects become not negligi- 
ble. This question has been the subject of extensive the- 
oretical activity over the past years, mostly concerned 
with the vibrational effects, for which we refer to a re- 
cent review^ and to the references therein. The role of 
the center-of-mass oscillations of a bridging site between 
the two leads has been well addressed, mainly via gener- 
alized Master equations in the context of nanoelectrome- 
chanical quantum-shuttle devices. On the contrary, the 
low-temperature quantum-coherent regime has been only 
slightly touched, and with rather controversial results at 
that. For instance, Al-Hassanieh et al^ made use of ex- 
act diagonalization procedure supplemented by a Dyson- 
equation embedding to conclude that conductance should 
be suppressed in resonance conditions for arbitrary cou- 
pling strength between the center-of-mass motion and 
the hybridization with the leads, and both at finite and 
vanishing charging energy. This result was questioned by 
Mravlje et ed^ who found, by a variational procedure and 
for finite charging energy, that the center-of-mass motion 
does not affect perfect transmission at resonance. 

In this paper we address the same class of questions, 
concerning the role of the center-of-mass motion at low 
temperature, in a different type of systems, namely 
metallic break junctions (BJ).** In a BJ the metal bridge 
or neck, initially forming a single solid body strongly 
bonded with the leads, is mechanically broken apart typ- 
ically at criogenic temperatures. The conductance drops 
prior to breaking typically takes place through a sequence 
of plateaus corresponding to thinning of the neck, down 
the ultimate monatomic contact, whose conductance is 



of the order of the conductance quantum Go = 2e? jh, 
where e and h are the electron charge and Planck's con- 
stant. These plateaus, are interpreted in terms of ballistic 
conductance, which in the adiabatic Landauer-Buettiker 
linear response thcory^'^ is controlled by the few residual 
one-electron conduction channels and by their respective 
transmittivity. The instant when the left and the right 
leads are separating, the physical bridge between the two 
is as a rule a single metal atom - as indicated by the last 
conductance plateau^ Here, non-adiabatic effects could 
in principle be caught right at the moment of separa- 
tion. The bridge atom, initially strongly attached to both 
leads, eventually detaches from one of them to remain af- 
ter separation exclusively attached to the other. In the 
process, the atom coordinate will move for a while in a 
double well effective potential. Therefore, between the 
initial solid metal-metal nanocontact, held together by a 
strong bond and with electrical properties governed by 
ordinary ballistic conductance, and the broken contact, 
there is room for a transient state where a new regime 
involving double well tunneling may be relevant. As the 
double well initially develops out of a flat single well, the 
two well minima can to a good accuracy be considered 
equivalent; moreover the barrier separating them is ini- 
tially very weak, which calls for quantum tunneling, even 
when the atom mass is not small. A skematic sketch of 
this system is shown in Fig. [TJ If the mechanical breaking 
takes place slowly enough in time, the dynamics of the 
bridge atom nucleus tunneling in the double well may 
be approximated by that of a two-level system (TLS), 
whereas the electronic level of the bridge atom, assumed 
to be nondegenerate, gives rise to a resonant electronic 
level. 

Within these assumptions, the physics might be as- 
similated to that of conduction electrons scattering off 
TLS's in bulk metals, proposed by Vladar and Zawad- 
owski - as a possible realization of a two channel Kondo 
(2CK) model.-^. This idea recurred several times in re- 
cent years in the context of a variety phenomena in met- 
0^^ ^9,10^11,12^13,14 ^ although again rather controversially. 
According to Aleiner et al^ in fact, the appropriate 
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FIG. 1: Sketch of a mechanical break junction 



high-energy cut-ofF of a TLS coupled to conduction elec- 
trons is not the electron bandwidth but rather the oscilla- 
tion frequency within each potential well, since above this 
energy the conduction electrons follow adiabatically the 
motion of the atom. Since the Kondo temperature, Tk^ 
is typically exponentially smaller than the high-energy 
cut-off, the conclusion of Aleiner et al. is that Tk is too 
small to be relevant. 

The BJ problem proposes now a new interesting phys- 
ical situation which we treat here in a slightly different 
model, arriving at interesting conclusions about the zero 
temperature conductance and its temperature evolution. 

By means of the Numerical Renormalization Group 
(NRG) r^ we calculate the influence of the two level sys- 
tem on the ballistic conductance across the bridge atom. 
We find that the zero temperature, zero voltage conduc- 
tance is always finite and fractional. However it is found 
to drop quite fast to zero as the temperature increases. 

The paper is organized as follows. In Section II we in- 
troduce the model Hamiltonian and discuss the param- 
eters chosen. In Section III we first solve some limiting 
cases by mean of analytical methods. In Section IV we 
study the low energy properties of our model by mean of 
NRG. We show that conformal field theory (CFT) pro- 
vides a strikingly direct interpretation of the low lying 
spectrum obtained by NRG. In Section V by using both 
our NRG routine and a Fermi liquid theory we give an 
estimate of the conductance of our model. In Section VI 
we finally summarize and comment our results. 



II. MODEL HAMILTONIAN 

The physics and language of our model is inspired by 
a bridge atom suspended between two one-dimensional 
metallic leads and moving quantum-mechanically in a 
symmetrical double well potential, although the model 
could equally other pseudospin variables coupled to a 
ballistic conductance channel. As a simplification we will 
assume that the dynamics of the atom nuclear coordinate 
is that of a TLS. We introduce a pseudo-spin variable Tz 
identifying the atom position, r^ = 1 and Tz = —I when 



the atom is in the minumum close to the right (R) and 
left (L) lead respectively. With this definition, assuming 
the bridge atom wavefunction to be real^, the quantum 
tunneling operator between the two wells corresponds to 
the Pauli matrix t^. 

The right and left metal leads are modeled as semi- 
infinite chains, site-label n — l,...,oo, with nearest 
neighbor hopping, amplitude —t, and creation (annihi- 
lation) operators cl^na i'^ana)^ whcrc a = R, L and the 
spin a =1 , |. The bridge atom is endowed with a single 
nondegenerate electronic orbital (the bridge level), of cre- 
ation and annihilation operators dj. and d^, respectively, 
constituting the ballistic conducting channel. The elec- 
tron hopping amplitude from the leads to the bridge level 
is assumed to depend explicitly on r^. When the atom 
is in the right well (r^ — +1), the level is more coupled 
to the R chain, amplitude — io(l -I- 7) with < 7 < 1, 
than to the L chain, amplitude — io(l — 7)7 and viceversa 
when the atom is in the left well (r^ = —1). Therefore 
the model Hamiltonian reads 



T^' — i /_^ / . 7 . ^anaCgn+la + H.C. 
Q=R,L o- n=l 

-toY. (1-7^.) [cl,^d^ + H.c. 

(T 

-A^r^-V^T:, ^ 4d^. (1) 



The last term represents the electron assisted tunnel- 
ing of the bridge atome nucleus arising from the influence 
of the atom's state of charge on the height of the barrier 
of the double well tunnelling potential^. In principle this 
type of assisted tunneling process includes other possible 
operators that couple the bridge level and the nuclear 
pseudospin coordinate, provided (given our assumption 
of a symmetric double well and equivalent leads) they are 
equally even under under reflection with respect to the 
center of the double well (we will call this even parity). 
The last term in ([T]) is therefore just one of the operators 
that presumably might possess a large matrix element, 
involving the bridge level charge occupancy. In later cal- 
culations below we will actually consider more general 
assisted tunnehng operators too. 

One can note at the outset that the model in ([T|) is 
closely related to a 2CK model, the role of the spin being 
played by the lead label, R and L, for the conduction 
electrons and by the pseudospin t that identifies the TLS, 
while the role of the silent channels is played by the real 
spin a. An alternative way of writing ([T]), which may be 
convenient in some cases, is by introducing the even (e) 
and odd (o) combinations 
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and formally defining 

through which the model ^ is rewritten as 
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are the local generators of the flavour SU(2), with ct" 
the Pauli matrices. 

Our model H in @ differs from the anisotropic 2CK 
model (O since 

• the even (fl") chain has one more site than the odd 
one (-IL); 

• in the even (ff) chain the hopping between sites 1 
and 2 differs from the others; 

• a local magnetic field Aj; acts on the pseudospin. 

In addition, 14 in (|4]) is generally coupled to an opera- 
tor more complicated than Tf , unlike J^ in ^. This 
difference has no effect when A^. = 0, in which case 
our model Eq. @ will display the conventional 2CK be- 
havior, but plays an important role when a finite A^; 
drives the model away from the 2CK fixed point. Specif- 
ically, we found that models with different ^ and and rj 
in Eq. (|3]), may fall into two different classes: 

(i) if 77 = and ^ = 1, hence the assisted tunneling 
term VxTx in ^ is proportional to Tf, or, more 
generally, if 



where 



Ve^V2to, Vo = V2to' 



(5) 



e = l + '7, 



(8) 



In iJ of (HI) we in fact included additional assisted tun- 
neling operators with coupling constants parametrized by 
T] and ^, which are missing in Eq. ([T]). In the even — odd 
formulation, the analogy with a 2CK model is much more 
explicit, especially once we rotate the pseudospin by 7r/2 
around the y-axis, even, e, and odd, o, labels playing the 
role of spin up (fl") and down (JJ-). A similar model was 
recently proposed by Zarand in the context of TLS's in 
metals jii according to whom the presence of the resonant 
level may push the equivalent 2CK model into a strong 
coupling regime with a large Kondo temperature of the 
same order as the high-energy cut-offl4. For comparison, 
we may also write the conventional two-channel flavour- 
Kondo model (after a 7r/2 rotation around the y-axis of 
the flavour pseudo-spin) 



then an intermediate 2CK crossover regime should 
survive in presence of a small but finite A^,; 

ii) if Eq. ([8]) is not satisfied, them this crossover regime 
is likely to be absent for any A^, 7^ 0. In this case 
the model with ^ = 77 = can be taken as repre- 
sentative of all the others. 

We note that the condition ^ means simply that the 
assisted tunneling operator, 

a 

is orthogonal to the local charge density. 



U-2CK = -^ XI X 51 ^"»'^^"»+i'^ + -^-^- 
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The fact that such a property discriminates betweeen 
two quite distinct classes of behaviors suggests that the 
charge degrees of freedom play in this problem an ac- 
tive role, unlike in conventional Kondo models, as we are 
going to discuss in what follows. 



III. PRELIMINARY ANALYSIS OF THE 
MODEL 

Simplifying the double well dynamics of the bridge 
atom to a TLS form permits a numerical analysis of the 
original model Hq H]). We performed that analysis by 
means of the numerical renormalization group ^^ and the 
results will be presented and discussed later. Prior to 
doing that we can, exploiting the analogy with a 2CK 
problem, discuss first some instructive limiting cases of 
([T]) that can be easily understood. 

First, if Aj; = V^ = the model describes a conven- 
tional electron hopping across the bridge level with in- 
equivalent leads because of 7 7^ 0. In particular, for any 
value of Tz , the zero temperature differential conductance 
in units of Go = 2e^/h is readily found to be^^ 



_G 

Go 



1 + 7' 



(9) 



If 7 = with finite A^, and T4, it is more convenient 
to use the even-odd representation in which the conduc- 
tance is 



_G 

Go 



sin^ (Se - So) , 



(10) 



where Se and So are the phase shifts at the chemical 
potential in the even and odd channels, respectively, de- 
termined by coupling of the leads to the bridge level. 
By solving the one-dimensional scattering problem and 
choosing for simplicity ry = ^ = 0, we find that 



A - ^ 1 !^ 
^ " 2 + 2<2 ^"' 



So - 0, 



Since the bare electron hopping to is of the order of the 
eV, which is many orders of magnitude larger than both 
Ax and Vx, one can safely treat the latter terms pertur- 
batively within the path-integral formalism originally de- 
veloped by Anderson and YuvaP'' and by Hamann^^ for 
the single-channel Kondo and Anderson-impurity mod- 
els respectively. That approach had in turn been built 
by extending the Nozieres-De Dominicis solution^'^ of 
the X-ray edge singularity to a succession of emission- 
absorption processes. In our problem, because of the 
presence of the silent spin-channel and of the bridge level, 
it is necessary to resort to a multichannel extension of the 
Anderson- Yuval formalism^^i^^, where the perturbation 
expansion consists of a series of pseudo-spin flips induced 
by the operator t^. What is important in the calculation 
is the phase-shift difference suffered by each channel at 
any pseudo-spin flip. In the present case the most conve- 
nient representation is in terms of R and L leads. We do 
not present details of the calculations, since as it turns 
out the final result can be inferred by very simple argu- 
ments. Because as was said when 7 = 1 and r^ = -1-1, 
only the R lead is hybridized with the level, while the 
L lead is untouched, R acquires a phase shift 5^ = 7r/2, 
corresponding to a resonant level model, while for the 



left lead, L, S]^ = 0. Viceversa, for r^ 

the L lead that is coupled hence S]^ 

Therefore the phase shift differences in the pseudospin 



-1, it is only 
while S2 = 7r/2. 



-1 to r^ 



1 are Sr 



7r/2 



flip from Tz = 
and Sl = Sj^ — S^ = —tt/2 for each spin a, which here 
plays the role of a silent channel. This is exactly the 
location of the so-called Emery-Kivelson pointf^^ which 
also coincides with the intermediate coupling fixed point 
of the 2CK modeP^'^*'. Interesting enough, this situation 
should also correspond to the maximum Kondo tempera- 
ture attainable,"'^* confirming Zarand's expectationil. We 
find that, at equilibrium, the perturbative expansion of 
the partition function coincides with that of the general- 
ized resonant level model 



so that 



G 44 



is always finite. 
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A. Asymptotic solution for strong electron-nucleus 
coupling: 7 = 1 

The parameter 7 in ([T]) measures the strength of 
"electron-phonon" coupling between the bridge atom and 
the leads. When the bridge atom double well is tiny, the 
two minima are close, and 7 will be small; in a wide dou- 
ble well, with the bridge atom very close to either R or 
L leads, 7 will be large (while A^ will correspondingly 
be small). The upper limit for 7 is 7 — 1, when the 
bridge atom in the left (right) well only couples to the 
left (right) lead. As it turns out, this limit is interesting 
by itself. 



7Y* = Ho [ipf.i'sf] - Ax 



^}(o)/ + /tv/o) 



V.\lY[f-f][^lfiO) + i',f{0)], (12) 

where F = 4i§/i is the hybridization width of the d- 
level, which plays the role of the high-energy cut-off, and 

Tio [■'/'/: 4'sf] is the continuum limit of a non- interacting 
Hamiltonian on a closed chain for two different chiral 
Fermi fields ipf{x) and ^Jsf{x) that move with Fermi ve- 
locity v,^^ namely 

'HoVPf,i^sf]=iv ^ dx-ipl{x)dx->Pai^)y 

a=f,sf 

and, finally, / and /t arc the annihilation and creation 
operators of an auxiliary fcrmion satisfying P f — 1/2 = 

Tz- 



Here we labeled the fields following Emery and Kivel- 
son^ to stress the fact that the role of spin s (here the 
real spin a) and of flavour / (here the R and L leads) 
are interchanged in our model with respect to the con- 
ventional 2CK model. 

Unlike the Emery-Kivelson Hamiltonian^^ for the 2CK 
model, in our case a pseudo-spin field A^; is present, 
which spoils the anomalous 2CK behavior.-'' For any fi- 
nite Ax , the spectrum of the Hamiltonian (|12p is Fermi- 
liquid like, corresponding in fact to a 2CK model in pres- 
ence of a magnetic field applied to the impurity, a case 
studied by Affleck, Ludwig, Pang and Cox.^^ We further 
note that the original Hamiltonian ([4]) is invariant under 
a generalized parity operator 



V=T, (-1) 



No 



(13) 



where No is the total number of electrons in the odd 
channel. Since a Fermi-liquid spectrum implies that the 
TLS - the Kondo impurity - asymptotically dissolves 
into the conduction bath, it follows that the value on 
each state of the generalized parity operator P^ turns 
effectively into the "Fermi- liquid" parity (—1)^°. This 
observation implies a zero-bias conductance dictated by 
the form (jlOp in the low-energy spectrum. 

We conclude by briefly discussing the other limit A^r — 
0, when the bridge atom is so heavy, or the barrier so 
large, that double well tunneling is suppressed. Here the 
model flows to the 2CK fixed point, and here it is well 
known that the elastic scattering ^-matrix at the chem- 
ical potential is zero.^* Since the even and odd channels 
correspond in our model to the spin up and down chan- 
nels in the 2CK problem, both have vanishing 5-matrix, 
hence the conductance is zero. For an infinitesimally 
small magnetic field acting on the impurity spin, it was 
shown in Ref. 23i that a Fermi liquid behavior is recovered 
with a phase shift difference of 7r/2 between the two spin 
channels. The translation of this result in our case is not 
so straightforward since, in the absence of any coupling 
to the TLS, i.e. Vo — Vx — 0, the even and odd phase 
shifts are flnite unlike the conventional 2CK. Actually, 
since the even chain has one more site, the "bare" phase 
shift difference is already tt/2. One possibility appears to 
be that the tt/2 phase-shift difference acquired by switch- 
ing on an inflnitesimal A^, at the 2CK fixed point adds to 
the "bare" value to give a total difference of modulous 
TT. This would imply zero conductance for A^, ^ Vx, rais- 
ing to nonzero by increasing A^,. We shall see that this 
is actually what happens if Eq. ([5]) is satisfied, namely if 
the assisted tunneling does not involve charge degrees of 
freedom. In the opposite case, the conductance behavior 
is more complicated. 



IV. NUMERICAL RESULTS 

We address the spectrum of the model Hamiltonian 
-ff of (HI) by standard NRG^^, whose results we are go- 



ing to present in this section. Tentative values of the 
Hamiltonian parameters which we adopted are a conduc- 
tion bandwidth 2Do ^ 2 eV, the attempt frequency D ^ 
10-2 eV,i^iii Ve ~ K ~ 0.2 eV and T4 ~ A^ ~ IQ-^ Ve- 
As discussed previously, only the conduction electrons 
with energy smaller than the attempt frequency are in- 
volved in the pseudo-spin screening. In order to enforce 
this condition, we take a flat conduction-electron density 
of states of bandwidth 2 eV (the chemical potential is 
zero), but we assume that only the conduction electrons 
with energy —D < e < D are coupled to the local degrees 
of freedom. Consequently, we perform the NRG proce- 
dure only on these electrons, which amounts to assume 
an effective bandwidth 2D r^ 2 x 10^^ eV, yet with a 
flat density of states equal to the original one, namely 
0.5 eV"^ = 0.5 X 10~^/D. Using the attempt frequency 
D as our energy unit, the net result in the Wilson chaiiJ^ 
is a renormalization of 




v: 



e(o) 



(o) 
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which keeps the d-level hybridization width at the chem- 
ical potential invariant, while Vx and A^; rescale trivially 
into themselves: 



Vx 

Vx^ —D, 



At 
D ' 



implying Vx — Ax ^ 10"^ K,— the values we assume 
throughout. Moreover, to better identify each state of the 
spectrum, in the numerical calculations we implemented 
the spin SU(2) symmetry, the charge U(l) symmetry and 
the discrete parity defined by Eq. p3|) . 

Following the discussion of Sec. II, we ran NRG calcu- 
lations for the two different implementations of the elec- 
tron charge assisted tunnelling, i.e. case (i), in which 
Eq. ([5]) holds with 77 = and ^ = 1, and case (ii) with 
77 = ^ = 0. 

In Fig[2]we show the NRG-flow for the Hamiltonian in 
EqlUfor both cases (i) and (ii) above with Ax = and 
with Ve = V; = 2 and Vx = 10"^ Ve- The energy of the 
lowest lying eigenvalues are plotted as a function of the 
number N of NRG iterations corresponding to an energy 
(temperature) scale ujn = DA~^^^ where A is the Wilson 
discretization parameter (we henceforth set A = 2). At 
large N, the spacing between the levels, their degeneracy 
and the disappearance of any difference bewteen even 
and odd iterations N (see for instance Ref. [271 ) is typical 
of a 2CK. These results are summarized in TablH and are 
consistent with the conformal field theory prediction^Sii^ 
for the 2CK. 

The numerical results clearly show that, whatever the 
form of electron assisted tunneling, the system has a 
2CK behavior at low temperatures. The Kondo temper- 
ature Tk is conventionally estimated as D A~^'^''~^''^, 
with Nc the NRG-iteration at which e.g. the first ex- 
cited state is 10% off its asymptotic valueJ^i^l We find 
that, while cases (i) and (ii) have roughly the same Nc, 



the latter is strongly influenced by Vo/Ve- In particular 
Vo/Ve — 1, namely 7 ~ 1, is an optimal choice that min- 
imizes Nc ^ 25, consistently with the previous analysis, 
and corresponds to a temperature of few hundredths of a 
Kelvin. Remarkably, even and odd iterations are hardly 
distinguishable after very few iterations. That seems to 
be a property of the 2CK model right at its fixed point 
- the fixed point with the highest Tk ~ Z?,^^'^^ - which 
would imply that the above estimate of Tk is a strong 
underestimation of the real one. However, we cannot ex- 
clude the possibility that the even-odd collapse of the en- 
ergy levels might simply indicate a preliminary crossover 
to a regime where the effects of Ve and Vo are fully es- 
tablished while those of Vx are still negligible. 
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FIG. 2: NRG-flow of lowest eigenvalues for the model in Eq.|31 
with Aa; = 0. The case (i) is analysed in the left panel, the 
case (ii) in the right one. 
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TABLE I: Lowest energy NRG spectrum of the Hamiltonian 
H of[4]for Aj, = 0. The energies E{NRG) are given in units of 
the fundamental level spacing and compared with the confor- 
mal field theory prediction E{CFT). For each eigenvalue we 
indicate its degeneracy (deg) together with its quantum num- 
bers Q (charge), S (spin) and parity P defined in Eq. (|13p . 



the quantum numbers of each eigenvalue differ substan- 
tially from that model. In the flavour 2CK model, label- 
ing states with Q, S and the flavour T (see Eq. ([7])), 
one expects the lowest energy spectrum of Table |lll 
This spectrum is determined within conformal field the- 
Qj, y30,3i ^y so-called fusion of the free-electron spectrum, 
to the left in Tab. |lll with the fiavour primary field with 
T = 1/2. 
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TABLE II: Lowest energy spectrum of the 2CK model © for 
Va ~ 0, left table, and Va 7^ 0, right table, as expected by 
conformal field theory. 

By contrast, we found that the NRG spectrum that we 
actually find. Tab. [H can be obtained starting from the 
2CK one in Tab.HIlin the following way. 

(1) First we decompose the fiavour SU(2) -^ U(l)xZ2, 
where U(l) stands for the free bosonic theory that 
represents the ^-component of the flavour fleld, and 
Z2 is an Ising conformal field theory (see §18.5 in 
Ref. I33) . This decomposition leads to the spectrum 
in Tab. 1111133 
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TABLE III: Lowest energy spectrum of the 2CK model upon 
decomposing the flavour SU(2) into U(l) XZ2. Tz is the quan- 
tum number that defines the U(l) theory, while Z2 corre- 
sponds to the coset theory, which is an Ising one. 



We note that, although the level spacings and degen- 
eracies are those of the conventional 2CK model Eq. © , 



(2) Next we shift the charge Q and z-component of the 
fiavour Tz by -|-li^ This corresponds to the fact 



that the even chain has one more site. In this way 
we obtain the spectrum in Tab. lIVl which coincides 
with that one in Tab. [H including the degeneracy 
of each eigenvalue. 
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TABLE IV: Lowest energy spectrum obtained from the one 
in Tab. Hn] upon shifting Q and T^ by +1. 



We note that, if we recombine the charge U(l) with 
the Ising to form an isospin (charge) SU(2) theory, the 
spectrum becomes equal to the conventional 2CK one in 
Tab. ni with the role of Q played by T^ and that of T 
played by the isospin. In other words, it seems that, al- 
though the original model is not invariant under isospin 
SU(2) symmetry, the fixed point does in fact recover that 
symmetry. This unexpected result is confirmed by the 
spectrum calculated during the renormalization group 
procedure. Indeed, after very few iterations, the ground 
state becomes and stays for all A^ > 1 doubly degenerate 
with quantum numbers {Q,S) = (+1,0), (—1,0). 

The above observation also clarifies why the charge 
degrees of freedom play an important role once Ax is 
turned on. As said, a finite A^^ is equivalent in the 2CK 
language to a magnetic field on the impurity site, which is 
known to be a relevant simmetry breaking perturbation 
destroying the anomalous 2CK behaviori^ Indeed, we 
find that, as soon as A^; 7^ 0, the spectrum flows to a 
Fermi-liquid one that can be interpreted as independent 
even and odd electron channels suffering different phase 
shifts 6e and Jo- 
in Fig. [3] we show the NRG flow of the low energy 
spectrum for A^; = 10~'*T4 for the two cases (i) and 
(ii). The asymptotic spectrum can be straightforwardly 
interpreted using the single-particle spectrum of Fig. [3] 
and combining all possible single-particle excitations. In 
particular, we find that, for very small A^ ^ Vx, Se — 
So = for case (i) and 6e — So — 7r/4 for case (ii). A 
difference between the two cases is apparent also in the 
way their approach to the asymptotic behavior. In fact, 
for the same values of A^; ^ T4 , a crossover region with 
a Ax = 0-spectrum is still visible in case (i) but not at all 
in case (ii), see Fig. [31 This different low energy behavior 
has its counterpart on the conductance behavior, as will 
be discussed in the next section. 

We conclude this part by emphasizing that for realistic 
A3: — V^ no crossover is visible in the spectrum, which 



might suggest the absence of any intermediate tempera- 
ture regime dominated by the singular behavior of the 
2CK fixed point. However, this statement should be 
taken with caution, since, as discussed above, the flow, 
even at Aa; = 0, is quite atypical and does not allow for 
a precise determination of Tk- Indeed, for Ax 7^ 0, it 
remains true that the spectra of even and odd iterations 
collapse very fast. However, unlike the case A^^ = 0, 
the levels at even and odd iterations with equal energy 
have opposite charge Q and parity P. This compares 
well with the role of a local magnetic field in the 2CK 
at its fixed point: levels at even and odd iterations with 
opposite spin quantum number Sz collapsed Therefore, 
although we tend to believe that the above estimate of 
Tk ~ 10~* D is correct, we cannot exclude that the ac- 
tual value could be much larger. 




FIG. 3: NRG flow of the lowest eigenvalues for the model 
in EqU] with Aa; = 10~^K;. Case (i) left panel, case (ii) 
right panel. Even and odd iterations correspond to solid and 
dashed lines, respectively. 



V. CONDUCTANCE 

We mentioned earlier that the zero-bias conductance in 
the 2CK state, A^; = 0, is zero because the scattering ma- 
trix of both the even and the odd channels are zero^^. For 
finite Ax , the recovery of Fermi- liquid behavior allows us 
to estimate the conductance by the difference Sg — So, see 
Eq. (|10p . which can be extracted by the spectrum, for 
instance by calculating the energy difference between the 
two lowest energy states with (Q, S, P) = (0, 1/2, 1) and 
(0, 1/2,-1) in units of the level spacing: 



Se- So ^ TT (E^ 



(0,1/24) ~ ^(0,1/2,-1)^ 



These two energies correspond to the cost of adding 
an even electron, {Q,S,P) = (0, 1/2, 1), or an odd one, 
{Q,S,P) = (0,1/2,-1), to the ground state, which has 
quantum numbers (—1, 0, 1), see Fig. H) 
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FIG. 4: Graphical representation of the Fermi-liquid spec- 
trum for Ax 7^ 0. The even(odd) single particle energy levels 
are equidistant, however the even spectrum is shifted with re- 
spect to the odd one. The ground state is obtained by filling 
each level below the chemical potential, _B = in the figure, 
and has quantum numbers {Q, S, P) — (—1, 0, +1). All possi- 
ble excitations can be generated by combining single-particle 
excitations. We show for instance the two lowest energy ex- 
citations that amounts to adding one electron, either even, 
(0,1/2,-1-1), or odd, (0,1/2,-1), which we use to evaluate 
the phase shift difference 5^ — So- 
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FIG. 5; Top panels. Conductance in units of the conductance 
quantum Go as a function of temperature for different values 
of A^/Vx ~ 1 for case (i) (left panel), and (ii) (right panel). 
Bottom panels: zero temperature conductance as a function 
of the ratio A^/Vx for model (i) (left panel) and (ii) (right 
panel) . 



We calculate this phase shift, hence the zero-bias con- 
ductance, as a function of the temperature T (extracted 
from the NRG iterations) , for different values of the ra- 
tio Ax/Vx- We note however that, while we are quite 



confident about the values at low temperatures, those at 
high temperatures must be taken with caution since the 
spectrum is still far from a Fermi liquid one. The results 
are shown in the top panels of Fig.[n]for case (i) and case 
(ii) and realistic values of A^/Vx ^ 1 (red-bold curves). 
In both cases there is a significant thermal crossover with 
very small conductance before the asymptotic low tem- 
perature regime is reached. At zero temperature, the 
conductance is zero if A^: = 0. However, as soon as 
an infinitesimal A^, is turned on, the zero temperature 
conductance stays in case (i) but jumps to Go/2 in 
case (ii), see bottom panels of FiglS] For realistic val- 
ues of Ax — Vx, the zero temperature conductance is in 
all cases finite, G ^ 0.5 -^ 0.9 Go, and smaller then the 
unitary value. 



VI. DISCUSSION AND CONCLUSIONS 

In summary, we have discussed the influence in the 
transport across a bridge atom of its quantum mechani- 
cal center-of-mass motion, whose dynamics in the double 
well case we have approximated by that of a two-level 
system-. In this regime, the two equilibrium positions 
of the bridge atom play the role of a pseudo-spin, whose 
dynamics is influenced by the electron hopping from the 
contacts into its valence orbital. This realizes effectively 
the same physics of a magnetic atom or a quantum dot 
bridging between two leads, the role of spin played by 
the position of the atom and the real spin playing the 
role of an additional flavour index. It is speculated that 
this hypothetical situation might be applicable to a metal 
break-junction caught right at the breaking point, when 
the central atom bridging the two contacts develops, al- 
though for a very short time interval, a double-well po- 
tential before collapsing finally onto one of the two. 

We find that, as long as the atom can tunnel between 
the two contacts, the zero-bias conductance at zero tem- 
perature is finite, although smaller than its value in the 
solid metal- metal nanocontact, with a single well for the 
bridge atom. This finite conductance seems at variance 
with the earlier result by Al-Hassanieh et al-, accord- 
ing to which the zero-temperature conductance at reso- 
nance should vanish at zero bias when the center-of-mass 
motion modulates the hopping amplitude into the leads. 
The discrepancy might be due to our two- level- system ap- 
proximation or, more likely, to the different low-energy 
accuracy of NRG with respect to the numerical technique 
employed by Al-Hassanieh et al-. Indeed, we find that 
the finite-temperature conductance, which should corre- 
spond to the effective zero temperature value obtained 
with less low-energy accuracy, decreases quite rapidly to- 
wards zero with increasing temperature. 

In the limiting (and unrealistic) case of a vanishing 
spontaneous tunneling, A^; = 0, in spite of a finite as- 
sisted one, Vx 7^ 0, the model displays a two-channel 
Kondo behavior, again with vanishing zero-temperature 
conductance. For finite Aa; <C Vx, the zero tempera- 



ture conductance is found either to remain zero or to 
jump to 1/2 of the unitary hmit (the conductance quan- 
tum), G = 0.5Go = e^//i, depending on the form of 
the assisted tunneUng. On the contrary, for reahstic 
vahies of A^; ~ Vx, the conductance is always finite, 
G- 0.5-^0.9 Go. 

A critical aspect of the model is that, with the realis- 
tic parameters used, distinct signatures of the two-level 
system dynamics could be hard to observe at tempera- 
tures around 4 K commonly used in metal break-junction 
experiments.^ Even harder could be the detection of 
possible manifestations of two-channel-Kondo anomalies. 
Cooling to lower temperature would offer the possibility 
to observe these effects. Time resolved conductance ex- 
periments could show the tunneling regime as a transient 
just before breaking and a coherent Kondo-like regime 
could be reached for light-mass shuttling-centers. For in- 



stance, hydrogen atoms or molecules moving onto and 
into mechanically controllable break junctions (see e.g. 
Refs. 34,35,36). In that case, the conductance plateaus 
found below the unitary limit, could be ascribable to the 
two-level system dynamics, similarly to that found in our 
model, see Fig. [SJ now shifted to higher temperature 
scales. A possible realization could be a metal contact 
bridged by a malone aldehyde molecule, where a hydro- 
gen bond is known to shuttle quantum mechanically be- 
tween two equivalent positions^!. 
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